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We develop an adiabatic Floquet picture in the length gauge to describe the dynamics of a 
hydrogen atom in an intense laser field. In this picture, we discuss the roles played by frequency and 
intensity in terms of adiabatic potentials and the couplings between them, which gives a physical and 
intuitive picture for quantum systems exposed to a laser field. For simplicity, analyze hydrogen and 
give the adiabatic potential curves as well as some physical quantities that can be readily calculated 
for the ground state. Both linearly and circularly polarized laser fields are discussed. 
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I. INTRODUCTION 

The dynamics of atoms in intense, ultrashort laser 
pulses involves the absorption and emission of many pho- 
tons, making non-perturbative methods necessary. By 
far the most common such method is the direct numeri- 
cal solution of the time-dependent Schrodinger equation. 
Unfortunately, this approach does not often provide the 
kind of simple physical picture that lets one understand 
the dynamics and predict what will happen under dif- 
ferent circumstances. Moreover, numerically solving the 
time-dependent Schrodinger equation accurately is, in 
most cases, a demanding task. These shortcomings have 
played a large role in promoting the three-step (or sim- 
pleman's) model [H, of ionization dynamics. In its nor- 
mal application, each step of the this simple model draws 
from a different theoretical formulation, although it can 
be rigorously derived with a few approximations Q. 

At the same time, the non-perturbative treatment of 
molecular dissociation in intense, ultrashort laser pulses 
has led to a different — but simple and useful - physical 
picture that guides many researchers' understanding of 
these systems. This picture combines the Floquet rep- 
resentation with the Born-Oppenheimer approximation, 
producing Born-Oppenheimer potentials that include the 
effects of the laser field. All of the insight earned in un- 
derstanding field-free molecular dynamics can thus be 
directly applied to molecules in a laser field. Impor- 
tantly, these "dressed potentials" provide a description 
of the laser-induced dynamics using only a single theo- 
retical formulation. These dressed potentials have also 
provided the interpretations behind the mechanisms of 
bond-softening, vibrational trapping (bond-hardening), 
above-threshold dissociation, and zero photon dissocia- 
tion. Recently, the picture has even been generalized 
to include ionization-dissociation channels, leading to 
the prediction of a new mechanism — above thresh- 
old Coulomb explosion — which was subsequently mea- 
sured 

A natural question to ask, then, is whether such a Flo- 
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quet approach can be applied to atoms in intense fields 
and, if so, does it yield a similarly useful physical pic- 
ture f?\. In this paper, we answer this question by gen- 
erating dressed potentials for the electron's motion in a 
hydrogen atom exposed to an intense laser. It is the gen- 
eration of dressed potentials that distinguishes our work 
from the substantial body of work using the Floquet ap- 
proach for atoms in intense fields 0. Our approach does, 
however, bear a closer resemblance to treatments of Ry- 
dberg atoms jt^l2i| and is very similar in spirit to the 
recent work of Miyagi and Someda 13;]. They empha- 
sized the comparison of circular and linear polarization 
in the velocity and acceleration gauges. While they use 
the adiabatic potentials for radial motion in the acceler- 
ation gauge to interpret their numerical results, we focus 
on developing the length gauge adiabatic potentials into 
a more general framework for understanding and pre- 
dicting the response of an atom to an intense laser pulse. 
We further test the quantitative predictive power of the 
potentials by calculating physical observables in the adi- 
abatic approximation. 



II. THEORETICAL BACKGROUND 

The Floquet theorem has been invoked many times in 
the past to solve the time-dependent Schrodinger equa- 
tion for an atom in an intense field. In fact, Floquet 
calculations provide some of the most accurate non- 
perturbative ionization rates available. With one pri- 
mary exception, though, no simple picture like that for 
molecules has emerged from this previous work. Floquet 
ideas did, however, provide an elegant understanding for 
that exception: the stabilization of atoms against ioniza- 
tion in a high-frequency laser field We will thus be 
applying the Floquet theorem in a non-standard way. 

For simplicity, and to allow the focus to be on the 
new representation we introduce, we will consider a CW 
laser field. This simplification has the additional benefit 
of allowing easy quantitative comparison with previous 
calculations. Our approach can be generalized to the 
case of a laser pulse using the ideas described in Q. 

The Hamiltonian for a hydrogen atom in an intense 
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laser field is given in atomic units by 



as 



H = 



rrV^ - - 



(1) 



df 



(5) 



using the dipole approximation and choosing the length 
gauge (see App. For a linearly polarized CW laser 
field, the electric field is 5 = focos(wt)z. In the above 
expressions, r and z are the electron's coordinates with 
the z axis chosen along the laser polarization; £o is the 
amplitude of the laser's electric field; and a; is the laser 
frequency. 

Since the Hamiltonian is periodic, H {1+2-11 / oj) = H{t), 
the Floquet theorem states that the solutions of the time- 
dependent Schrodinger equation 



have the form 



H'ip{r,t) 



(2) 



(3) 



where (j>{r,t) = (j>{r,t + 2tt/uj). Upon substitution into 
Eq. ([2]), one finds the eigenvalue equation 
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This equation shows that, in a sense, the Floquet ap- 
proach allows us to treat time like a coordinate. For 
later convenience, we define the effective Hamiltonian % 



The primary reason that the molecular Floquet ap- 
proach described above is so useful is that it produces 
easy-to-interpret potentials that include the effects of the 
laser field. So, our goal is to similarly produce dressed 
potentials for an atom. The most natural way to do this 
is to take r to be an adiabatic parameter and solve 

Had*. = f/.(r)$.. (6) 
In this case, the adiabatic effective Hamiltonian is 



1 



d 



■Had = TT^ — 



2r 



dt 



(7) 



where is the squared orbital angular momentum op- 
erator. The eigenstates from Eq. ([5]) form a complete set 
at every r so that the wave function cf) can be written 
without approximation as 



E 



1 



(8) 



The eigenvalues Uy{r) of Eq. (|6]) are precisely the po- 
tentials we seek. That they serve this purpose can be 
seen by substituting Eq. ([8|) into Eq. (|4]) which yields the 
following equations for F,j{r) 



J 



1^ 

2(ir2 



Mr) - 



2 ^ 



dr 



F,,{r)^eF,{r). 



(9) 



These equations are identical to those one obtains in 
the Born-Oppenheimer representation, and the deriva- 
tive coupling matrices P and Q are defined as usual by 



P..'(r) = ((<i>.|-|<i>.' 

g..'M-((<i>.i^l<i>. 



(10) 



The double bracket notation indicates that these matrix 
elements should be integrated over one cycle of the field 
in time as well as over the electron's angles. Note that the 
coupled radial equations ^ , including the potentials and 
coupling, are time-independent. And, in the absence of 
the field, the equations decouple, yielding the usual sep- 
arable solutions for a hydrogen atom. It is worth noting 



that this adiabatic representation is, in principle, exact. 

To solve the adiabatic equation (|6]) , we expand the spa- 
tial degrees of freedom in $ on the spherical harmonics. 
We can guarantee the time periodicity of $ — and thus 
oi (j) — by expanding its time dependence via the discrete 
Fourier transform. We thus have: 



(11) 



nlm 



explicitly showing the parametric dependence on r. The 
Floquet index n ranges from — cx) to oo; Z is the angular 
momentum; and m is the projection of the angular mo- 
mentum on the z axis. In this basis, the matrix elements 
of T-Ladir) for a linearly polarized laser field are 
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{{nlm\'Ha.d\n'l'm')) = 6nn'Su'6r, 



1(1 + 1) 



2r2 



nUJ + {Sn,n' + 1 + Sn,n'-i)-^r£Q{Yi„i\cos9\Yi,jn')- (12) 



The term {Yim\ cos 9\Yii m') can be evaluated analyti- 
cally in terms of 3j symbols. The diagonal terms 
((nZm|'Had|"-^"i)) are the diabatic potentials. 

We can simplify the potential curves by using the 
dipole selection rules implicit in Eq. (IT2|) . Taking the ini- 
tial state to be the field-free hydrogen ground state and 
choosing it to correlate with n = in the field, only the 
diabatic channels that couple to the (n, Z,m) = (0,0,0) 
channel are relevant. Therefore only those diabatic chan- 
nels with both n and I even or n and I odd are allowed. 
Further, all channels with non-zero m can be eliminated, 
so this label will be suppressed. The physical meaning 
of the Floquet index is now clear, it represents the net 
number of photons absorbed or emitted by the hydrogen 
atom in the laser. 

Before proceeding to the results and discussion, we 
note that by selecting the length gauge the whole matter- 
field interaction is included in "Had since it is local in 
r. In the velocity gauge, however, this is not the case 
since the matter-field interaction contains a term with 
a derivative in r which is not included in Had- While 
Had in the acceleration gauge does include all the terms 
involving the field, it is also not ideal for reasons that 
are discussed in App. [Xj To be clear, even in the length 
and acceleration gauges in which all terms involving the 
field are included in T^ad, it is not true that all effects of 
the field are included in the adiabatic potentials — the 
non-adiabatic couplings are also field-dependent. From 
this discussion, it is clear that the adiabatic potentials 
are gauge-dependent, as are the non-adiabatic couplings. 
However, this gauge-dependence is just a matter of the 
distribution of physical information between potentials 
and couplings; any physical observable calculated in any 
gauge will be properly gauge-independent if calculated 
exactly, i.e. by including all of the potentials and their 
couplings. 



III. RESULTS AND DISCUSSION 
A. Linearly Polarized Fields 

1. Potential Curves 

Figure [Ija) shows the diabatic potential curves, in- 
cluding only the lowest two Is for each n. The diabatic 
potentials are thus partially dressed in that they have 
been shifted by the net number of photons exchanged 
with the field. They do not, however, reflect any further 
effects of the field as all such effects lie in the couplings in 
this representation. Passing to the adiabatic representa- 




FIG. 1; (color online.) (a) The diabatic Floquet potentials 
for a hydrogen atom at uu — 0.2 a.u. in a linearly polarized 
laser field. The thick blue curve denotes the initial channel, 
(b) The corresponding adiabatic Floquet potentials at, £o = 
0.05 a.u. The Floquet blocks included here are —10 < n < 10 
and the angular momentum states included are 1 — 0,1, 2, 3. 



tion by diagonalizing Had produces fully dressed poten- 
tials that include much more of the effects of the laser. 
These adiabatic potentials are shown in Fig. [TJb), and 
each crossing of diabatic potentials from Fig. [IJa) has 
now become an avoided crossing. A transition between 
diabatic channels represents the absorption or emission 
of a photon or photons (according to An) as does pass- 
ing through an avoided crossing on the same adiabatic 
channel. Transitions between adiabatic channels, which 
occur predominantly at avoided crossings, may or may 
not involve the exchange of photons with the field. 

Although Fig. [IJb) has more potentials than the typi- 
cal molecular example, the parallels between the atomic 
and molecular cases are rather clear. For instance, tak- 
ing the initial state to be the field- free hydrogen Is 
state (indicated in both panels of Fig. [T] at an energy 
of —0.5 a.u.), we see from Fig. [IJb) that the electron can 
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tunnel through the barrier induced by the field to emerge 
on potentials correlating to the n=2 and 3 thresholds. In 
the process of tunneling, the wavepacket thus absorbed 
three w=0.2 a.u. photons. 

If the tunneled wavepacket reaches r — ^ cxd on one of 
the n—3 channels, then it has ionized — the analog of 
bond-softening in molecules. The final kinetic energy of 
such an electron is the difference between the initial en- 
ergy and the final threshold energy, which in this case 
is 0.1 a.u. Since the electron absorbed two photons 
to ionize, dipole selection rules dictate that it should 
end up with I— I or 3. /^From the diabatic potentials 
in Fig. [ija), we can identify the most likely pathway as 
being (n,0 = (0,0) (1,1) (2, 2) '-^ (3,3) by 
following the crossings. Comparison with the adiabatic 
potentials in Fig. [ijb) show that these are the channels 
that form the barrier through which the electron tunnels 
to ionize. We thus predict that the ionized electron will 
primarily have / = 5. Physically, then, we expect the 
angular distribution of photoelectrons to have primarily 
/-wave character. 

If the wavepacket ends up on the n—2 channel instead, 
then the electron has simply been excited since the n—2 
threshold energy lies above the initial energy. To get to 
the n=2 threshold, the wavepacket must emit a photon 
at the crossing around r—5 a.u. Based on the small gap 
at that crossing compared to the kinetic energy at that 
r, the electron will most likely traverse the crossing dia- 
batically and thus ionize. Should it emit a photon and 
remain bound however, arguments similar to the ioniza- 
tion case predict that the excitation pathway should lead 
mainly to 1=0 rather than 1 = 2. 

These adiabatic potentials thus provide considerable 
qualitative — and quantitative — information. In fact, 
within this picture, we can also distinguish simultane- 
ous absorption or emission processes from sequential pro- 
cesses. Along the excitation pathway described above, 
for instance, the process is initiated by the simultane- 
ous absorption of three photons during the tunneling, 
followed after some delay by the emission of one pho- 
ton. "Simultaneous" in this case means that the photons 
were exchanged with the field at a single avoided cross- 
ing. The three-photon absorption during tunneling really 
takes places over an r range of about 1.5 a.u. as the po- 
tential has clear n=0 character only up to about r=2 a.u., 
and 71=3 character only after about r=3.5 a.u. The elec- 
tron, of course, takes a finite time to cross this distance. 
This example is rather extreme since most avoided cross- 
ings are much more localized as can be seen in Figs. [T] 
andm 

Being able to identify where the transitions occur in 
r could be useful for the attosecond control of the elec- 
tron's dynamics. The travel time of the electron between 
crossings can be estimated, for instance, making it pos- 
sible to identify opportunities for control either through 
the pulse length or through the delays between multiple 
pulses. Such control based on the Floquet potentials has 
been demonstrated in 4]. Control is also possible. 



in principle, for processes having two or more pathways 
to the same final state. In this case, some mechanism 
to control the phase accumulated along each pathway is 
needed. 

FigureOshows the adiabatic potential curves for differ- 
ent laser frequencies and field strengths. These potentials 
show field strengths corresponding to 3.5x10^^ W/cm^ 
(£"0=0.01 a.u. [15]), 8.8x10^3 W/cm^ (£o=0.05 a.u.), and 
3.5x10^"* W/cm^ {£q=0.1 a.u.). The range of frequencies 
represented, w=0.2-0.6 a.u., mean that the number of 
photons required for ionization range from three to one, 
respectively. In our adiabatic Floquet potentials, then, 
ionization of II(ls) proceeds either by tunneling through 
a barrier or by passing over it. 

In addition to barriers, the adiabatic potentials in 
Fig.[2]also have wells. The well most likely to play an im- 
portant role in ionization of II(l,s) forms the other half of 
the avoided crossing that produces the barrier discussed 
above. Its minimum lies between r=l and 2 a.u., depend- 
ing on uj, but is energetically accessible at -0.5 a.u. only 
for a;=0.4 and 0.6 a.u. in the figure. When it is energet- 
ically accessible, there is a chance that part of the initial 
wavepacket will be trapped in the well — the analog of 
molecular vibrational trapping. Should this happen, it 
might be observable as a reduced ionization rate. Or, 
in a time-dependent measurement, it might be observ- 
able as ionization on different time scales as the trapped 
portion of the wavepacket should take longer to ionize. 
Other analogs to molecular phenomena can be identified, 
some of which will be described below. 

The discussion of these adiabatic Floquet potentials 
would not be complete without some comment on their 
range of applicability. Even though the adiabatic Flo- 
quet representation described here and implemented in 
Eq. ([9]) is exact, the potentials themselves form a use- 
ful qualitative picture only when there are not too many 
of them. Physically, this condition is clearly achieved if 
the adiabatic parameter r is indeed the "slowest" vari- 
able in the system. Practically, so long as the photon 
energy is not too much smaller than the ionization en- 
ergy, then the adiabatic picture will be useful. For II(ls), 
the a;=0.2 a.u. case shown in Figs. [2lJa-l)-(a-3) is about 
as small as is useful. Decreasing uj further not only re- 
quires a larger range of n, but also a larger range of I 
since Z > |n| follows from the dipole selection rule. The 
complexity of the picture can thus grow rapidly. 



2. AC Stark Shifts 

The adiabatic Floquet potentials can also be used to 
make nontrivial quantitative predictions without full- 
blown calculations. One such example is the AC Stark 
shift. Strictly speaking, in a laser field, all bound states 
become resonances — a point made clear within the 
present representation. The shift in the position of that 
resonance from the field-free value is the AC Stark shift. 
Since our adiabatic potentials include the effect of the 
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FIG. 2: The adiabatic Floquet potentials for a hydrogen atom in a hnearly polarized laser field at (a-1) u = 0.2 a.u., 5o = 
0.01 a.u.; (a-2) lj = 0.2 a.u., fo = 0.05 a.u.; (a-3) w = 0.2 a.u., fo = 0.05 a.u.; (b-1) w = 0.4 a.u., £^ = 0.01 a.u.; (b-2) 
— 0.4 a.u., ~ 0.05 a.u.; (b-3) oi = 0.4 a.u., — 0.1 a.u.; (c-1) oj = 0.6 a.u., = 0.01 a.u. (c-2) w = 0.6 a.u., fo = 0.05 a.u.; 
(c-3) Lk) = 0.6 a.u., i^o = 0.1 a.u. Seven Floquet blocks, n = —3 to 3, and four angular momentum states are used for these 
curves. 



laser field, we will estimate the ground state energy in 
our initial channel using a sirnple WKB formalism (in- 
cluding the Langer correction (l6|). 



/' 



dr\ 



1 

2 ' 



(13) 



which is applicable only when the shifted state is still 
below the barrier shown in Fig. [T] Here, ri and r2 
are the classical turning points, and eis is the shifted 
ground state energy. Because the initial adiabatic po- 
tential in Fig. [1] is, in fact, quite complicated with very 
sharp avoided crossings at r values smaller than the bar- 
rier position, we used a diabatized potential V{r) in the 
calculation (see Fig. [3]). This approximation can be phys- 
ically justified by the fact that these sharp crossings are 
most likely traversed diabatically. 

The ground energy shifts as a function of field strength 
for several lu are shown in Fig. |4l For comparison, we 
also include the results calculated by the non-Hermitian 
Floquet matrix method 17[ , which predicts a larger shift 
than ours in all cases. Surprisingly, our results give nearly 
perfect quadratic behavior of the shifts with £q in the 
regime where the barrier is still present. Quadratic be- 
havior is expected for the shift from second-order pertur- 




FIG. 3: (color online.) The diabatized potential curve (thick 
dashed blue line) supporting the AC Stark-shifted ground 
state. Here, lu = 0.2 a.u. and Eo = 0.05 a.u. 



bation theory. The differences between our results and 
the non-Hermitian Floquet calculations, which should be 
quite accurate, are most likely due to our neglect of the 
non-adiabatic couplings in Eq. (HI 
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Eo(a.u.) 
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FIG. 4: The absolute value of ground state energy shifts of a 
hydrogen atom in a linearly polarized laser field as a function 
of field strength for (a) At oj = 0.18 a.u. (b) At uj ^ 0.19 a.u. 
(c) At [J = 0.20 a.u. and (d) At ui = 0.22 a.u. Note the actual 
shifts are negative. The dotted line shows Sg for comparison. 
The vertical line indicates the field at which over-the-barrier 
ionization begins. 



3. Ionization Rates 



is lower than the ground state energy when ui > 0.25, ui 
needs to be larger than 0.25 a.u when the system is ini- 
tially in the ground state, for the transition formula to 
be applied. 

Strictly speaking, to get the ionization rates we need to 
include all the adiabatic potentials and the couplings be- 
tween them, which would be tough numerically because 
of the sharp avoided crossings. Fortunately, we can treat 
most of the sharp crossings as purely diabatic. We also 
simplify the potential curves by seeking the most impor- 
tant pathways relevant to the physical process we are 
interested in. Here we select sequential one-photon pro- 
cesses, hence neglecting simultaneous multi-photon ab- 
sorptions (emissions) . 

Since Landau-Zener only provides the the transmission 
probabilities, the flux of electrons through the crossing is 
needed to calculate the ionization rates. We determine 
it from the classical argument that an electron with 1=0 
and an energy of —0.5 a.u. passes the crossing point twice 
when it moves from one end of its trajectory to the other. 
The frequency of the electron passing through the cross- 
ing point is thus Q = l/7ra.u. The ionization rate is then 
obtained as following: The rate of ionization to the n-th 
Floquet channel is then 




n=0 



■n=l 



-11=3 



r (a.u.) 

FIG. 5: (color online.) The pathways for a ground state 
wavepacket in a linearly polarized laser field. The horizon- 
tal arrow indicates the initial state at an energy of —0.5 
a.u. (dashed line). Here u — 0.4 a.u. and £o = 0.05 a.u. 

Another quantity that can be estimated from our adi- 
abatic potentials is the ionization rate. The ionization 
rates are calculated for the cases where the ground state 
energy is higher than the top of the barrier shown in the 
above adiabatic potential curves. We use the Landau- 
Zener transition formula for curve crossings to ob- 
tain the transition probabilities. The frequency regions 
where our calculations of ionization rates could apply are 
restricted by the prerequisites of the Landau-Zener for- 
mula: the energy of the wave-packet is high above the 
curve crossing point and only two-curve crossings are in- 
volved. Figure [5] shows the corresponding pathways for 
the present purpose. Since the first crossing point (be- 
tween diabatic channels Z = 0,n = and / = l,n = 1) 



n n ^'.'+1 (1 - ^' 



(14) 



vJ=0 



taking into accounr the absorption of each of the n pho- 
tons. The Landau-Zener transition probability is 



Ti,i+i 



1 — exp — 



27r 



1/2 



(15) 



V;_/+i is the diabatic coupling, -Fi.;+i is the classical 
"force", and is the "velocity". They are all eval- 

uated at the crossing /+i = y^T+l)/a; and are given 
by 



V, 



1,1+1 



■l + l 



3 LJ 



{Yio\Yio\Yi^i^o)So, 



F, 



_2iJ_ 

VT+T' 




(16) 



The ionization rates for each threshold shown in the adi- 
abatic potential curves for uj = 0.6 a.u. in Fig. [5] are 
calculated by Eq. The outgoing kinetic energy of 

the ionized electron is simply obtained by finding the dif- 
ference between the initial energy (ground-state energy 
of a field-free hydrogen atom) and the outgoing threshold 
—nuj. 

Figure [B] shows the ionization rates for multi-photon 
absorption up to n = 10. Since in our picture the laser 
field is purely monochromatic, the outgoing kinetic en- 
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ergies can only have the discrete values En = nuj — eis, 
which corresponds to the peaks of the above threshold 
ionization (ATI) spectrum. It is seen from Fig. [6] that 
the rates drop nearly exponentially as n increases, which 
is the same feature seen in ATI experiments [19]. We also 
converted our total ionization rates to the half-width of 
the shifted ground-state energy for a rough comparison 
with the data from [l3|. It can be seen in Table|T]that our 
results are on the same order of magnitude as the results 
from Ref. [l7j , which is very good agreement considering 
the simple scheme we used. 



FIG. 6: The ionization rates as function of photoelectron en- 
ergy at u = 0.6 a.u., £o = 0.1 a.u. for a linearly polarized 
field. The dashed line indicates a pure exponential decrease 
for comparison. 



TABLE I: Comparison of ionization rates for H(ls) at a; = 
0.6 a.u. in a linearly polarized field. 



£o(a.u.) 


F/2'^(a.u.) 


F/2''(a.u.) 


V2-0.01 


0.3937 X 10-5 


0.1253 X 10-3 


72-0.025 


0.1414 X 10~3 


0.7836 X 10-3 


V2-0.05 


0.5618 X 10-3 


0.3144 X 10-2 


V2- 0.075 


0.1250 X 10-2 


0.7107 X 10-2 



"From present calculation 

''From non-Hermitian Floquet method [17| 



B. Circularly Polarized Field 



Our approach can be generalized to circularly polarized 
laser fields. In this case, we take z axis of the electronic 
coordinate to be perpendicular to the laser polarization, 
such that S = £{)[cos{Lut)x -I- sin(wi)y]. The matrix ele- 
ments of Had {f) are now 



{{nlm\'Had\n'l'm')) = Snn'Sii'S„ 



l{l + l) 1 

^ — nuj 

2r^ r 



m I 



r So- 



il?) 



Again the number of potentials curves we need to include 
can be simplified by selection rules. As with the linearly 
polarized laser filed, only those diabatic channels with 
both n and / even or n and I odd are allowed. For cir- 
cularly polarized laser field, however, m is now restricted 
by m = n. 

We show the diabatic and adiabatic potentials for a 
circularly polarized laser field in Fig. [T] There is a signif- 
icant difference in these the potentials from the linearly 
polarized case: at any value of r, there is always a low- 
est potential curve no matter how large n becomes. This 
is certainly not true in for linear polarization since in- 
creasing n always adds a lower energy potential. The 
difference stems from the fact that in a circularly polar- 
ized field m — n which makes the angular momentum 
I increase as n increases. Since the diabatic potentials 
are quadratic in I but linear in n, they become increas- 



ingly repulsive for larger n. As illustrated in Fig. [Tl^b), 
we observe that at large distances, the lowest potentials 
become parallel to each other and diverge to — oo like 
—r^. The electron will therefore accelerate to infinity 
if it stays in one adiabatic channel. We interpret this 
unphysical situation by failure of the aiabatic approxi- 
mation for large r in the circularly polarized field. We 
first note that in our picture, the change of the adiabatic 
potentials signifies the gain of electronic energy from the 
laser field in the region of large r where the Coulomb po- 
tential is negligible. This energy gain can be understood 
by the following simple classical argument: the way we 
calculate the adiabatic potential curves is equivalent to 
confining the electron to a fixed r and finding its station- 
ary motion. In a circularly polarized field, the stationary 
motion is just that the electron moves in a circle in con- 
stant speed with the electric force always perpendicular 
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r (a.u.) 

FIG. 7: (color online.) (a) The diabatic Floquet poten- 
tials for a hydrogen atom in a circularly polarized field with 
oj — 0.2 a.u. The thick blue curve supports the initial 
state, (b) The corresponding adiabatic Floquet potentials 
for So — 0.05 a.u. The Floquet blocks included here are 
— ll<n<ll and the angular momentum states included 
are < Z < 11. 




r (a.u.) 

FIG. 8: (color online.) Rescattering dynamics for (a) linearly 
and (b) circularly polarized light. The red arrows indicate 
the ionization pathway most relevant to rescattering, while 
the blue dashed curves trace the corresponding potential. In 

(a) , |n| < 15 Floquet blocks were included with Z < 5; in 

(b) , |n| < 21 and I < 21. In both figures, uj = 0.2 a.u. and 
£o = 0.05 a.u. 



to the direction of the motion. The eletron energy for 
large r is then given by in agreement with the 

leading order of |?7(?')|. The physical classical motion of 
the electron during ionization, however, is a spiral tra- 
jectory where the asymptotic energy of the electron does 
not change drastically. Therefore, to correctly describe 
the motion of the electron in the adiabatic Floquet pic- 
ture, a lot of adiabatic channels need to be coupled for 
large r, which makes pathway analysis impractical. 

In much the same way as for linear polarization, we 
can calculate the AC Stark shift of the H(ls) ground 
state in a circularly polarized field. Now though, the 
lowest adiabatic potential can be directly used without 
diabatization. In Fig. [HI we show the AC Stark shifts 
calculated by Eq. (IT51) . Similar to the linearly polarized 
case, the shifts follow an £q scaling. Comparing with 
Fig. mja) for the linearly polarized case, we see that our 
approximation predicts the AC Stark shifts for circular 
polarization are larger than for linear polarization, even 
when compared at equal intensities, i.e. Sun = V^Sdrc- 

When one-photon ionization is dominant, we can cal- 
culate the ionization rates in a similar manner as dis- 
cussed for linearly polarized fields in Sec. IIIC. The po- 
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FIG. 9: (color online.) The absolute value of ground state 
energy shifts of hydrogen atom in circularly polarized laser 
field as function of field strength at a; = 0.18 a.u.. The dashed 
line shows £0 ^ comparison of the behavior of the shifts, the 
vertical line shows the field where above of barrier ionization 
occurs. 
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tentials and transitions are identical to those but with dif- 
ferent couphngs. The expression for the ionization rates 
is then still given by Eqs. (fTSj) and (fT6|) . with Vi^i+i re- 
placed by 



Vi,i+i = ^Jy^-^(Yi+U+i\Yii\Yu)£o- (18) 

The first diabatic coupling Vq i has the same value as in 
a linearly polarized field if the field So is scaled by 1/ v^, 
which correspond to the field with the same intensity. 
Since the total ionization rate for linearly and circularly 
polarized fields in the perturbative limit are dominated 
by the single-photon process, they should be related by 
a simple rescaling of the field strength. 

Having potentials for both linearly and circularly po- 
larized light, one interesting question to consider is 
whether they can clearly describe the difference in rescat- 
tering behavior between the polarizations. Namely, elec- 
trons ionized in a linearly polarized field can be driven 
to rescatter from their parent ion; electrons ionized by a 
circularly polarized field do not rescatter. To answer this 
question, we show in Fig. 8 the adiabatic Floquet poten- 
tials for both polarizations. For linear case in Fig. 8(a), 
the electron tunnels through the barrier to ionize. Some 
part of this ionizing wavepacket can find itself on the 
n = 2, ? = curve after a few transitions. This potential 
is attractive, and this part of the the wavepacket gets 
reflected back towards the ion as expected. For the cir- 
cular case in Fig. 8(b), however, once the electron passes 
over the barrier to ionize, it encounters only repulsive 
potentials and thus never rescatters. 

IV. SUMMARY 

In this paper we have developed the adiabatic Floquet 
picture to describe the dynamics for hydrogen atom in 
monochromatic laser field, although the analysis is triv- 
ially extended to any effective one-electron atom. The 
radial distance of the electron r is adopted as adiabatic 
parameter, leading to a set of potential curves that in- 
clude the effect of the laser field. The potentials are com- 
pletely analogous to Born-Oppenheimer potentials can be 
interpreted in the same way. In particular, very simple 
approximations such as WKB and Landau-Zener can be 
applied to make semi-quantitative predictions from these 
potentials that are non-perturbative in the field strength. 

Just how fruitful this picture is remains to be seen. 
The diabatic curves have the distinct advantage of being 
very easy to generate, and already allow the application 
of simple approximations. While we did not explore it 



here, one of the more intriguing possibilities is using these 
curves to estimate the times at which particular photon 
transitions occur. Few other approaches provide this in- 
formation, especially in such a simple manner. Having 
established the reliability of our picture in this paper, we 
hope that such timing information can prove useful in 
controlling attosecond dynamics. 



Appendix A: Gauge issues 

While physical observables must be gauge invariant (in 
the limit of no approximations), the adiabatic Floquet 
potential curves are not. The choice of gauge is usu- 
ally made for computational convenience and clarity in 
analysis. There are, in principle, an infinite number of 
choices, but in practice three main gauges are usually 
discussed: (i) the velocity gauge, which is the minimal- 
coupling Schrodinger equation in the dipole approxima- 
tion; (a) the length gauge, which was the starting point 
for this paper; and (Hi) the so-called "acceleration ga.uge" 
(also known as the Kramers- Henneberger frame [20|,|2l|). 
As mentioned above, the velocity gauge is not attractive 
for an adiabatic representation because not all terms in 
the Hamiltonian involving the field are contained in "Had- 
The length and acceleration gauges, however, are both 
good candidates for an adiabatic picture since the laser 
field is completely included in Had (see also [l^)- Here 
we will briefly discuss the adiabatic Floquet picture for 
atomic hydrogen in the acceleration gauge. 

The adiabatic Hamiltonian for hydrogen in the accel- 
eration gauge is 

L'^ 1 d 

"Had = - ] 77TT~*^' (^1) 

2r"' |r — a(t)| ot 

where a{t) is the classical trajectory of a free electron in 
an oscillating field. In the present case of a CW laser, 

a{t) ~ zao cos{ujt) (A2) 

with QfQ = So/uj'^, assuming q;(0) = ci(0) = 0. 

One advantage of the acceleration gauge is that the 
effect of the electric field vanishes in the limit r ^ oo. 
Understandably, then, the usual field-free solutions are 
recovered, whereas in the length gauge the laser-atom 
interaction diverges linearly in r. A disadvantage of the 
acceleration gauge is the difficulty presented by the os- 
cillating Coulomb singularity. 

Applying the same adiabatic Floquet analysis as in 
Sec. ini we obtain the following matrix elements for the 
adiabatic Hamiltonian: 
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\nkrri\'Hi^d\n'k'm')) = Snni / sm6u1,{9) 

Jo 



sm 



6'u*,(e')e-*("-"')"*Ufc(6') 



Jo \/r^ — 2raQ cos{ujt) co8{9) + cos(wt) 



dtd9, 



r 
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FIG. 10: The adiabatic Floquet potentials for atomic hydro- 
gen in the acceleration gauge. The crosses indicate the ini- 
tial curve the system starts on when the electric field is zero 
(ao = 0). 19 Floquet blocks were used in this calculation 
(centered around n = 0). 



At that electric field strength and frequency, ao = 1.25, 
and in the figure the complicated nature of the curves 
is quite clear for r < 1.25. While in the acceleration 
gauge, the channels decouple at large r, they are strongly 
coupled at small r. Further, unlike the length gauge, the 
familiar dipole selection rules are no longer valid, and in 
fact all of field-free angular momentum states are coupled 
together near r = ap. This makes tracing a path through 
the potential curves very difficult, even locating the curve 
the system initially sits on is a challenge. The crosses in 
Fig. [TU] trace out the initial state when ao = 0. In the 
underlying mess of curves, one can see the semblance 
of the avoided crossings mentioned in Fig. [IJb), but the 
hopes of finding a path that passes though these crossings 
is daunting. While the thresholds are easy to determine 
in the acceleration gauge, the steep trade-off is in the loss 
of physical intuition in the region close io r — a^. 



where n and n' are again the Floquet indices, and Uk{9) 
replace the spherical harmonics as the basis functions for 
the 9 dependence. In our calculations, we took Uk{9) to 
be b-splines [22l|. 

In Fig. [TUl wc plot the adiabatic Floquet potentials in 
the acceleration gauge, using the same parameters as in 
Fig. mb), namely w=0.2 a.u. and £o=0.05 a.u. 
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